Functional macroscopic behavior of weighted ran- 
dom ball model 



O 

(N 

(N 

P? 
CM 

-I— > 



> 
in 
oo 

ON 

(N 
O 
O 



Jean-Christophe Breton and Clement Dombry 

Laboratoire MIA, Universite de La Rochelle, 17042 La Rochelle Cedex, France 
E-mail address: jcbretonOuniv-lr.fr 

Laboratoire LMA, Universite de Poitiers, Teleport 2, BP 30179, F-86962 Futuroscope- 
Chasseneuil cedex, France 

E-mail address: clement.dombryOmath.univ-poitiers.fr 



Abstract. 

by 



We consider a generalization of the weighted random ball model defined 



M(y) = / 



y -x 



N(dx, dr, dm) 



where N is a random Poisson measure on M. d x K + x R with a product heavy 
tailed intensity measure and h : K d — > K is a fading function . This functional can 



serve as a basic model for transmission with fading effect. In lBreton and Dombry 



(|2009h . he authors proved the convergence of the finite-dimensional distributions 
of related generalized random fields under various scalings and in the particular 
case when h is the indicator function of the unit ball in R d . In this paper, tightness 
and functional convergence are investigated. Using suitable moment estimates, we 
prove functional convergences for some parametric classes of configurations under 
the so-called large ball scaling and intermediate ball scaling. Convergence in the 
space of distributions is also discussed. 



1. Introduction 

We consider weighted random balls in 
N\ on W l xl+xl with intensity 



l d generated by a Poisson random measure 



n\(dx,dr,dm) = \dx.F{dr)G{dm) 

where A e R + and F, G are probability measures on K + and K respectively. For 
each 3-uplet (x, r, m), x represents the center of the Euclidean ball 5(x, r) and r 
its radius, m stands for the weight of the ball. The parameter A is interpreted as 
the intensity of the balls in M d . Su ch models ar e used for instance to r eprese nt a 



spatial communication network, see iKai I (|2006l ). lYang and Petropulu I (|2003l ). In 



this case, x represents a statio n transmitting a signa l, r the range of emission and m 



the intens ity of the signal, s e elBreton and Dombry 

particular lKai et al. I (|2007l) , IKai and Taaau I (|2007|) . Following this interpretation, 



(2009) and reference therein in 
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the signal m transmitted by x is received in some y S M. d if and only if y € B(x, r), 
and the quantity of signal received from the stations at y is given by 

/ ml B ( Xir )(y)N\(dx,dr,dm) = / ml B (o,i)((y-x)/r)N\(dx, dr, dm). 

(1.1) 

From a modeling point of view, it is natural to consider that the signal transmitted 
by x and received in y fades when y gets away from the station x. In order to take 
into account this phenomenon, we introduce a fading function h replacing 1_b(o,i) 
in (jl.ljl . more precisely the faded signal received at y from x is m/i((y — x)/r) and, 
assuming moreover that no interference occurs between the stations, the quantity 
of signal received at y is now given by 



M(y) = / mh((y -x)/r)N x (dx,dr,dm). (1.2) 

JR^xR+xR 

From a physical point of view, it is natural to assume that h(y) is a radially non- 
increasing function with h(0) = 1, < h(y) < 1 and limi|y|i_> +00 h(y) = 0. The 
function h is said to be radially non increasing if for all y S K d , the function 
r h4 h{ry) is non-increasing on [0,+oo). However, on a mathematical point of 
view, more general assumptions will be enough, see assumption QA3D below. 



In the sequel, we are more generally interested in the contribution 

M(ji) = f M(yHdy) 



of the model in a configuration of points y represented by a measure /x. For instance, 
the configuration reduced to the point y is represented by \x = S y and in this case 
M(5 y ) = M(y). It is natural to consider finite positive measures /1 on R d but 
our study supports signed measures fi with finite total variation. In the sequel, 
we shall note M. the set of such measures and we recall that, equipped with the 
total variation norm ||/x||jVf = ^t|(K d ), Ai is a Banach space. Actually i n order to 
make our study easier and in contrast with iBreton and Dombry" (2009), we shall 



consider signed measures [i with density, i.e. n(dy) — 4>(y)dy for some G i x ( 
Setting T Xir h(y) = h((y - x)/r) and (i[f] = j Rd f{y)fi(dy) for / e L 1 ^,^), the 
Fubini theorem allows to rewrite 

M{p) = I m^[T Xir h]N x (dx,dr,dm). (1.3) 

jR d xR+xR 

Note that the stochastic integral in (|1.3j) is well defined and the change in the order 
of integrals is authorized when 

/ \m/i[T x , r h] I n\(dx, dr, dm) 

il d xl+ xR 

< A|^|(M d ) (j \h(x)\dx\ (^J r d F(dr) \ (J \m\G{dm) \ < +00. (1.4) 

We will always suppose that all three integrals in (jl.4l) above are finite (precise 
assumptions on F, G and h are given in the set of conditions (A) below). Further- 
more, in this case, the expected value of M(p) is given by 

E[M(n)] = A^(K d ) I [ h(y)dy) ( f r d F(dr)) ( f mG(dr 
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In order to investigate the macroscopic behavior of the generalized random field 
(M(/z)) jU g.M, we apply a scaling x H> px, with p < 1. The scaling contracts the 
space and is interpreted as zoom-out in the model. Note that when p > 1, 
th e scaling becomes a dilatation of R rf and is interp reted as zoom-in. In contrast 



tn e sc aling becomes a dilatation ot IK ana is interp reted as zoom-m. in contrast 
to Bierm et al. ( 2010f) . Breton and Dombrv ( 20091 ). we focus in this note only on 



zoom-out (see below for further comments o n the relation between this note and our 
previous paper iBreton and Dombrv I (l2009h l. In order to derive non-trivial asymp- 



totics, the intensity A of the Poisson measure is adapted to the scaling procedure 



by allo wing A := X(p) to depend on the zooming factor p, see lBreton and Dombrv 
(2009). !^ote that the natural intensity X(p) corresponding to the scaling x H> px 
is X(p) = p~ d \. After the scaling, the generalized field becomes 



Mp(fi) = m(x[T x>r h]N Pi x(p)(dx.,dr,dm) (1.5) 

Jl d xl+xl 

where N p \( p } is the Poisson random measure with intensity 

n p,\(p)(dx.,dr,dm) = \(p)dx.F p {dr)G{dm) 

and F p is the image measure of F under r H> pr. We are finally led to investigate, for 
a proper normalization n(p), the fluctuations of the rescaled and centered random 
field 

M p (p) = n{p)-\M p { l j) - E[M p (p)]). 

In order to derive non-trivial asymptotics for the model (|1.5p . the distributions F 
and G driving the behavior of the radius r and of the weights m, and the shape 
function h must satisfy some conditions, denoted by conditions (A), that we state 
now precisely: 

• The probability G is assumed to belong to the normal domain of attraction 
of the a-stable distribution S a (a,b,r) with a £ (1,2], i.e. if X\, . . . ,X n 
ar e i.i.d. accord ing to G, n~ 1 / a {Xi + • • • + X n ) => S a (a, b, r). According 



to lFeller 1 11966, XVII. 5), this is equivalent to the following estimate on the 



characteristic function tpc of G: 

<p G (9) = l + i6T-<T a \6\ a {l + ibe{8)t&n(ira/2)) + o(\8\ a ) as 6 -> 0. (Ai) 

In the case a € (1, 2), a typical choice for G is an heavy-tailed distribution 
while for a — 2, G may be any distribution with finite variance. Observe 
that since a > 1 , G has a finite moment of order 1 . 

• The probability F is assumed to have a regularly varying tail satisfying 

r+oo 

F(r) := / F(du) ~ r ^>+oc Cpr^ 13 for some d < (3 < ad. (A 2 ) 

J r 

Here and in the sequel, f(r) ^r^+oo g(r) indicates that lim r _,. +00 /(r) / g{r) = 
1. Observe that, under ((A2J, the expectation of the volume of the random 
balls J R+ r d F{dr) is finite (see Lemma IA.3I below) and the bound (|1.4[) 
indeed holds true. 

• The shape function h is assumed to be continuous almost everywhere and 
such that 

h*(x) :=sup{|h(rx)| : r > 1} e L 1 {R d ) D L a (M. d ) . (A 3 ) 
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Note that this implies that h € L 1 (R d )DL a (R d ) and that if h is radially non- 
increasing, then h* = h. Indeed, h* is the smallest radially non-increasing 
function dominating h. 

The convergences of t he finite-dimensional d istributions (fdd) of M p were (essen- 
tially) already derived in Breton and Dombry ( 20091 ) under three different regimes 



depending on the behavior of X(p)p /3 (the so-called large, small and intermediate 
ball regimes). In this note, we actually focus on the corresponding functional con- 
vergence for the generalized random fields I M p (/i) ) . Since there is no natural 

functional space in which the random function /j, S M. i-> M p (fi) belongs (at least 
heuristically), we choose to consider a special parametric sub- family (//t)teK*' °f 

Ai and to investigate the tightness of the random fields ( Af„(«t) ) ■ Our main 

V /tea" 

contribution is to prove tightness in the space of continuous functions C(K P ) of such 
random fields under suitable conditions on the family (pt)tGVLP- As a by-product, we 
also obtain weak convergence of the random fields M p in the space of distributions. 

Let us conclude this intro duction with some comments on the relations between 
this note, our p r eviou s paper [Breton and Dombrvl (|2009[) and other related pa pers 
iMikosch et all (|2002t ). lKaj and Taqqu I (|2007l ). In lBreton and Dombry I (|2009h . fdd 



convergences are obtained for (M p (fj,)) fl when p belongs to some special subspace 
M. a ,p on which 

\p(B(x,r))\ a dx < C(r p Ar q ) for some p < /3 < q (1.6) 

d 

(roughly speaking the condition requires a control of the measures u (B(y:,r)) of 



both l arge and small balls, uniform in the centers of the balls). Moreover . iBreton and Dombry 



( 20091 ) deals simultaneously with the macroscopic behavior (i.e. p — > and F has a 



power law behavior in +oo of order j3 > d) and microscopic behavior (i.e. p — > +oo 
and F has a power law behavior in of order f3 < d). In comparison, in this note 
we deal only with the macroscopic behavior (i.e. p —> and f3 > d) and for special 
measures p{dy) = (j>{y)dy G M with density (j> e L 1 (R d ) n L a (R d ). But, F is 
not any more assumed to have a density and we only assume the tail condition 
QA2D . Moreover a shape function h is added in the model (| 1 . 2[) to take into account 
the fading of the signal in the communication network. But our main contribution 
in this setti ng is to derive tightne s s to o btain functional counterparts of the con- 
vergences in IBreton and Dombry (2009). In the particular case of the dimension 
d = 1, the model is related to the inf inite Poisson model in [ Mikosch et al. I ( 2002 ) 
and to the continuous flow reward in lKai and Taqqu I ( 20071 ). In both papers, the 



authors deal with the asymptotic behaviour of the random process corresponding 
in our setting to (M p (l[ ^(y)dy))t>o and the issue of tightness is adressed in both 
papers. 

The rest of the note is organized as follows: the main results are stated in 
Section [2] and proved in Section [3l Technical results are postponed in an Appendix. 

2. Main results 

First, we recall the finite-dimensional convergence for the generalized random 
field M p = n(p)^ 1 (M p (p) — K[M p (p)]). Actually, we state a slight modification of 
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the main results in lBreton and Dombrv 1 fl2009h replacing 1_b(o,i) therein by a shape 



function h satisfying QA 3 p . We shall abusively write p G L 1 (M d ) n L a (R d ), instead 
of p(dy) = (fi(y)dy with (fi G L 1 ^) n L Q (R d ). The symbol ^4 stands for the /dd 

y 

convergences and the symbol is also used throughout to indicate a functional 
convergence in the functional space X (for instance, in Theorem 12. 2\ X = C(M. d ), 
the space of continuous function of R d ). 

Proposition 2.1. Suppose conditions (A) hold. 

(1) (Large ball regime) If X(p)p@ — > +oo, then, setting n(p) — (A(p)p^) 1 / Q , we 
have as p — > 0: 

M p (p) M Z a (p), p G L 1 (R d ) n L Q (R d ) 
where Z Q is the stable field 

Z a {p)= \ p[T x , r h]M a (dr,dx) 

with respect to the a-stable measure M a with control measure a a Cpr~ Y ~^ drdx 
and constant skewness function b, where cr and b are related to G by QAi| ) . 

(2) (Intermediate ball regime) If X(p)p l3 — > a for some a g (0, +oo), then, 
setting n(p) = 1, we have as p — > 0: 

M P (M) ^ Ja(/x), /i G i 1 ^") n L"(R d ) 
where J a is the compensated Poisson integral 



Ja(n) = / m^[T x , r /i]A^ /3ia (dx, dr,dm) 

with respect to the compensated Poisson random measure Np a with inten- 
sity aCpr~ l3 ^ 1 dx.drG(dm). 

In the sequel, the finite-dimensional results are strengthened into functional con- 
vergence for a parametric sub-family of measures pt{dy) = </>t(y)dy,t g R p , in 
L 1 (R d ) n L Q (R d ), so that the generalized random field M p induces a p-dimensional 
random field (M p (/x t ))teR»- To that aim, we investigate the tightness in C(R P ) of 
this induced p-dimensional random field. The proof relies on a Censov criterion 
and on moment estimates for increments presented in Section 13.21 Actually in our 
setting, the relevant increments are generalized increments on blocks defined as fol- 
lows. Let s, t G W be such that s% < U, 1 < i < p, and consider the corresponding 
block [s,t] = n?=iN 3 *i]- The dimension of [s, t] is given by the number of indices 
i such that s% < ti. The generalized increment of a random field X = (X t )tsRj> on 
a block [s, t] in W is defined (up to a factor ±1) by 

X{[s,t\) :=^(-l) p -£^ e *X( Sl + ex{t x - Sl ),- ■ • , s p + e p (t p - s p )) (2.1) 

e 

where the sum above runs over e G {0, 1} P such that et — if Sj = ti. Similarly for 
measure pt with density (fit, we define the increment on a block [s,t] of (fit by 

</>[s,t] ■=^2(~l) P ^= iei <l>s 1 +e 1 (t 1 -s 1 ),...,s p +e p (t p -s p ) (2-2) 
e 

where again the sum runs over e G {0, 1} P such that q = if Sj = U. Such general- 
ized increments are easy to handle in our context and Example 12.21 below confirms, 
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for uniform-type measures, that such increments make sense in our setting. In order 
to control such increments, we introduce our main condition on the densities <pt of 
/it: we say that the family of densities (</>t)te»> satisfies property (P 7 ) for 7 > 1 if 
for all T > 0, there exists some constant Ct > such that for any [s, t] C [— T, T] p , 

11? <ct n i**- a *i- w 

i:Si <ti 

Here, ||</>|| 7 stands for the L 7 (R p )-norm of </>. The following examples justify that 
condition (P 7 ) is natural. 

Example 2.1. Let (/^t)teK* be the family of (signed) uniform measures on the blocks 
I '*] = rii=i [°) *«] : more precisely, <j> t = sign(ii) • • • sign(f p )l[ ,t]- Then, we verify 
that for any non-degenerated block [s,t] C R p , 0[ s ,t] — l[s,t] almost everywhere so 
that 

w<t>[s,t]\\!= n 1** — 

l<i<p 

In the case of a degenerated block [s, t] C [— T, T) p , we have 

= n 1**1 n \ti- s i\< c T n i**-**i 

i:Si—ti i:Si<ti i:Si<ti 

with Ct — max(l,T) p . Hence (P 7 ) holds true for all 7 > 1. Such uniform-type 
measures are u s ed to analyze cumulatiye wor kload in one-dimensional model, see 
Mikosch et aT~l (|2002l ). iKai and Taaau I (|2007t ) . 



Example 2.2. Let (fj,t)tem.p be a family of measures with densities <f>t such that for 
all I C {l,...,p} and all T > 0, 

sup <9/</> t (y) < +°° ( 2 - 3 ) 



te[-T,T]p 



7 



where 9/ is the differential operator defined for / = , . , , ik} by dj — d k /dti 1 ■ ■ ■ dti 
The following condition (P') (that implies condition (P 7 )) is satisfied: for any 
T > 0, there exists some constant Ct > such that for any [s,t] c [— T, T] p , 

11^11} <c T n iti-siV- (K) 



This fact is justified in Section 

We now state our results for the functional convergences in the large ball and 
intermediate ball regime. We recall that the limits Z a and J below are defined in 
Proposition 12.11 



Theorem 2.2. Suppose conditions (A) hold. Let (^(dy) — (f)t(y)dy,t <E K p , be a 
parametric family of measures in i 1 (R d ) n L a (M. d ) satisfying conditions (Pi) and 

{Pa). 

(1) (Large ball regime) If X(p)p 13 — > +00, then, setting n(p) — (A(p)p /3 ) 1 / Q , we 
have as p — > 0: 
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(2) (Intermediate ball regime) If X(p)p@ — > a > 0, then, setting n(p) = 1, we 
have as p — > 0: 

M^t)^J«W, ter. 

As a by-product of the moment estimates used to prove tightness, we obtain 
Holder-regularity properties in the case a — 2. This is the content of the following 
result. 

Proposition 2.3. Suppose conditions (A) hold with a = 2. 

(1) If G has a finite variance and the family of measures pt(dy) — (j)t{y)dy 
satisfies (Pi) and (P2), then the Gaussian limit process (Z2(pt))teR" 01 
Theorem E21 is 7-Holder for all 7 < 

(2) If there is k > 2p such that h £ L k (M. d ), G has finite moment of order k, and 
the family of measures pt{dy) — 4>t(y)dy satisfies (Pi) and (Pfc), then the 
limit process ( J a {pt))temp of Theorem l2.2l is 7-Holdcr for all 7 < 3( ^~' 3 — |. 

Observe that in dimension d = 1, we recover at the limit the fractional Brownian 



motion obtained at the limit in lKai and Taqqu I (2007) with the Holder-regularity 
7 <(3-/3)/2 6(i,l). 

Finally, we consider functional convergence in the space of distribution. Let T>(R d ) 
be the space of smooth compactly supported functions, and V(R d ) be its dual, i.e. 
the space of distributions. We show that, for all p > 0, M(p) can be seen as a 
random distribution and state functional convergence in T>'(R d ). 

Theorem 2.4. Suppose conditions (A) hold. 

(1) For each p > 0, M p induces a random distribution, i.e. the linear form 

V(R d ) -» M 



M p 



M p {4>(y)dy) 



is almost surely continuous. 

(2) (Large ball regime) If X(p)p 13 — > +00, then, setting n(p) — (X(p)p /3 ) 1 ^ a in 
M p (p), we have as p —> 0: 

M p {p)ZE^ Za {p), peV(R d ). 

(3) (Intermediate ball regime) If X(p)p^ — > a > 0, then, setting n(p) = 1 in 
Mp(p), we have: 

M» ^3 J ( M ), M GP(M d ). 

Remark 2.5. The result is still true if £>(R d ) is replaced by C^-(R d ) the space of 
compactly supported functions of class C k on M. d and V(M. d ) is replaced by the dual 
ofC k K (R d ), keN\{0}. 

3. Proofs 

The proof of Theorem 12.21 as usual for functional convergences, consists of two 
arguments: fdd convergences and ti ghtness. The first one is giv en in Proposition 
12.11 ; this is a slight modification of iBreton and Dombrv I (|2009l) and the changes 



are discussed in Section 13.11 The proof of tightness is given in Section 13.31 and it 
relies on moment estimates previously obtained in Section 13.21 Holder regularity 
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also relies on moments and cumulants estimates and Proposition 12.31 is proved in 
Section 13.41 Section 13.51 is devoted to the proof of functional convergence in the 
space of distributions V'(R d ). 

Recall that throughout the paper, we assume that conditions (A) hold, and in 
particular we consider d < (5 < ad. All the asymptotics are considered as p — > 0. 



3.1. fdd convergences. The results of Breton and Dombry ( 20091) do not apply di- 
rectly since the model we consider here is not exactly the same: in lBreton and Dombry 
(|2009t ). the probability F is assumed to have a density, the fading function h is re- 
placed by Ib(o.i) an d the fdd convergence is proved on the space M. a ,p of measure 
p satisfying (|1.6p . Here the radius distribution F and the shape function h are 
more general (see conditions (A)), but we obtain convergence only on the space of 
measures /tx(rfy) = 4>{y)dy with (f> S L 1 (R d )nL a (R d ) . However , a s brie fly explained 
below, verbat im changes in the proofs of Breton and Dombrvl ( 2009h shows that 
the results of iBreton and Dombry | (|2009t ) still apply in our context. 

Indeed, first a careful reading of the proofs in lBreton and Dombry I fj2009h shows 
that the existence of the density / of F is used only in Lemma 3 .2 therein. But this 



(2007) deriving the same 



lemma can be replaced by Lemmas 2 and 3 inlKai et al. . , _ „ 

result as Lemma 3.2 in lBreton and Dombry ( 20091 ) under the weaker assumption 
( IA2I ). Observe in particular that F„(l) = F(l/p) ~ C rp 13 and that the continuity 



requirement in Lemmas 2 and 3 in fkai et al. I (|2007| ) is ensured in our setting by 



condition ((A3]) and Lemma IA.2I below. 

Second, the bound (|1.6|) can be replaced by the following condition (|3.1[) : for 

h G L 1 (M. d ) n L a (R d ) and fi(dy) = 4>{y)dy with G L 1 (R d ) n L a {R d ), we have 



\lt[T Xtr h]\ a dx < C{r d Ar ad ), 



(3.1) 



this is justifies in Lemma I A. II in the Appendix with 7 := a therein. 

As a consequence, Theorems 2.4 and 2.11 in lBreton and Dombry ( 20091 ) remain 
true in our setting and ensure the fdd convergences stated in Proposition 12. II 



3.2. Moment estimates. As we will see, our results on tightness and Holder regu- 
larity strongly rely on moment estimates for the rescaled random field M p . Since 
the following properties of the moment are also interesting in their own right, they 
are stated in the following proposition. 

Proposition 3.1. Suppose conditions (A) hold. 

(1) Let < 7 < a. There exists some constant C :— C(F,G,h,a, f3,-f,d) 
depending only on F, G, h, a, /3, 7, d (i.e. not on p and </>) such that for all 
H(dy) = cj>(y)dy G i 1 (M d ) n L a (R d ) and all p > 0, we have 

■y/a 

E 



M p {n) 



< C 



<p) a 



(3.2) 



2 and G has a finite moment of order k 6 N \ {0, 1}, 
d ). There exists some constant C := C(F, G,h,a, (3,k,d) 



(2) Suppose that a 
h e L 1 (R d )r\L k 

depending only on F,G,h,a, /3,k,d such that for all p(dy) = <fi(y)dy £ 
L 1 ^) n L k (R d ) and all p > 0, 

Xp)p 



Ck(M p (p)) 



< C- 



n(py 



k(f3-d) k(kd-{3) 
I (fc-l)d || ^j| (k-l)d 

\k 



(3.3) 



Weighted random ball model 



9 



where Ck(M p (fi)j is the cumulant of order k of M p (p). 

As a by-product of these moment estimates and the finite-dimensional conver- 
gence (Proposition I2.ip , we obtain the following result stating the convergence of 
moments: 

Corollary 3.2. Suppose conditions (A) hold and p(dy) = 4>{y)dy with <f> G L 1 (R d )n 
L a (R d ). Let < 7 < a. 

(1) If \{p)pP — > +oo, then, setting n(p) = {\{p)p l3 ) 1 / a in M p (p), we have as 
p^O: 



E 



M p ( M ) 



E[\Z a {n)\ 



(2) If X(p)p 13 — > a, then, setting — 1 in M p (/i), we have as p — > 0: 



E 



m p (ao -^E[|j ( M )r 



In the case a = 2, if furthermore G has a finite moment of order ri > 2 and 
L 1 ^) n L n (R d ), then the above convergence of moments holds for all 
< 7 < n. 



Proof of Proposition [3711 
First point: The estimate for E 



relies on the following expression of the 



fractional moment fsee lvon Bahr and Esseen (1965) or lKai and Taqqu I (|2007l Eg. 
(60))): if A is a random variable with characteristic function fx{t) = E[exp(iiA)], 
then we have, for 1 < 7 < 2, 



+00 



E[|An=A( 7 ) / (l-l^^l 2 )^ 1 -^ 



(3.4) 



where 



A(rf) = (J (1 - cos^aT 1-7 *^ < +c 



Since M p (p) is a Poisson integral, its caracteristic function is given by (see Lemma 

USD 

(n)W = eX P ( / *G (n(p) -1 ^[Tx 1 r'i]) X(p)dxF p (dr] 

with * g (m) = / E (e mm - 1 - ium)G(dm). Hence, 

i-l^wWI 2 



1 — exp 



2Re(* G (n{p)- 1 9p[T^ r h]))\(p)d X F p (dr) 



< 1 - exp -2 



R d xR+ 



tf G (n(p)-^M[r x , r /i])| \{p)d*F p (dr 
2C(G)n(p)- a \(p)\6\ a [ W^ r h]\ a d*F p (dr) 



< 1 — exp ^— 

Using Lemma lA. 41 in the Appendix with 7 := a > /3/d, we have 

\»[T x , r h]\ a dxF p (dr) </£(</>) 



(3.5) 



(3.6) 
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where C(<f>) is given by 

and C(F) is a constant depending only on F. Plugging the bounds (I3.5[) and p.6[) 
in (13.41) . we obtain 

E 



M P ( M )| ] < A( 7 ) j (l-exp(-2C(G)n(p)-«A( /9 )/|0| Q C(0)))0- 1 -^0 
= A{i)A{a,i)(^pj'' {2C{G)C{4>)V' a (3.7) 



with a straightforward change of variables in (|3. 71) and 



A(a,7) 



/■+oo 

/ (1 - exp (-0°)) _1_7 d0 < +00. 
Jo 



This gives the result (|3.2[) with the constant 
C(F,G,h,a,P,~/,d) 



(a-l)0d 



(ad-p)a \ 7/« 



The case 1 < 7 < a is proved. The case < 7 < 1 comes from the previous 
case applied to any 7' G (l,a) combined with the Jensen inequality which implies 
E[|X| 7 ] < E[|X| 7 '] 7 / 7 '. 

Second point: We use the general form of the cumulant of a Poisson integral as 
recalled in Lemma [A. 51 Using Lemma [A.4I with 7 :— k > 2 > /3/d, we have: 



Hp) 

n(p) 



[ m k G(dm)) ( [ [ f i[T^ r h]} k dxF p (dr)) 
Jr J \Jm d xR+ ' ) 



5 ^^ g!t-m-^ m ^ mmi ^ ' lm ' a(im) 

This gives the bound ()3.3j) with the constant 
C(F,G,h,a,/3,k,d) 

(h — ~\\Ril Hfi-d) (kd-H)k / r- 

= V-w-o """ 1 ' 1 "' 

□ 

Proof of Corollary [3T2l 

We use the following basic result: if X n weakly converges to X as n — > +00 and 



E 



\x„p 



is bounded, then the convergence of moments E [|JT„| 7 ] — > E [|X| 7 ] holds 

for any < 7 < 7' (in this case, the family \X n \ J is indeed bounded in L 1 1 1 and 
hence equi-integrable) . As a consequence of this result, the finite-dimensional con- 
vergence stated in Proposition 12.11 and the moment estimates obtained in Proposi- 
tion l3.1l vield the convergence of moments for < 7 < a and < 7 < n respectively. 
The case 7 = n should be proved separately but we omit the details. □ 
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3.3. Tightness. This section is devoted to the second step in the proof of Theorem 
12.21 i.e. tightness. More precisely, we show that: 



Proposition 3.3. Under the assumptions of Theorem 
(M p (/i t )) teR P,p < 1, is tight in C(W). 



the family of random fields 



The proof of Proposition 13.31 relies on a suitable control of the moment of the 
generali zed increments of M a (p,t) wi th the follow ing Censov criterion f o r whi ch we 
refer to iDeshaves and Picard I ([1984 p. 16) or to iBickel and Wichura I (|197ll ) . 



Proposition 3.4. Let (X„)„ e n be a sequence of random fields on W such that: 

(1) The family of random variable (X„(0)) n£ N is tight; 

(2) For all T > 0, there are constants 7 > 0, 6 > 1 and Ct > such that, for 
all [s,t] c [-T,T]p and n € N, 

E[|X„([s,t])n <C T 1] \U~ Sl \ s . 

R emark 3.5. Proposition 13.41 i s a weaker version of the original criterion given 
in IDeshaves and Picard (1984) that we have adapted to our setting. Note that 
rii-Si/ti l*< — s i\ is the fc-dimensi onal Lebesgue measure of [s, t ], where k is the true 
dimension of the block [s,t]. In IDeshaves and Picard I (|1984l) . the Lebesgue mea- 
sure is replaced by general Radon measure with diffuse marginals. Originally, the 
criterion is expressed in terms of bounds on the tails P(|X „,([s,t1)| > x) instead 
of the moments E [|X„([s, t])| 7 ]. Moreover, it is explained in Deshaves and Picard I 
(|1984l ) that when fc-dimensional faces are involved, it is enough to check the condi- 
tion on blocks [s,t] such that Sj = tj — for the degenerated dimensions. For the 
sake of clearness, we do not insist on such general conditions. 



Proof of Proposition 13.31 

We apply Proposition 13.41 to the family X p (i) = M p (pt), p < 1. To that aim, 
observe first that the tightness of X p (0) is a consequence of the one-dimensional 
convergence stated in Proposition Furthermore, using the definitions (|2.1j) and 
(|2.2p together with the linearity of the generalized random field M p , we see easily 
that 

-MM]) = M p(A*[ s ,t]) with M[s,t](dy) = <P[ s ,t](y)dy. 

Proposition 13.11 -1) gives the following moment estimate for < 7 < a: 



E[|X p ([s,t])P 



= E 
C 



n(p) a 

Using (Pi) and (P a ) for <^[ s ,t]j we obtain: 



E[|x p ([ S ,t])n < c 



liP-d) 

(Q-l)d 



ll0M]||a? Il0[s,t]||i' 



-l(ad-/3) 



c 





j/a ^ 


n(p) a 




\Kp)p ' 


7/a . 


n{p)<* _ 





iW-d) , T(ad-g) 
Q ( Q -l)d" 1 " (a-l)d 




(3.8) 



(3.9) 
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Now, observe that the sequence \(p) p@ / n(p) a is bounded since under the two 
asymptotics investigated \(p)pP /n(p) a converge to some finite constant as p — > 0. 
Furthermore since d < (3 < ad and 1 + a ~ j3/d> 1, the exponent — (1 + a — 0/d) 
is (strictly) larger than 1 for 7 close enough to a. This proves the second condition 
in Proposition I3~4l Proposition 13.31 and thus Theorem 12.21 easily follow. □ 

We finish this section with the proof of Example 12.21 where the differentiability 
condition (I2.3[) is stated to be sufficient for (P^). 

Proof for Example 12.21 Let s < t. We first assume that, for all 1 < i < p, 

-T < si < U < T. We have 



0M] (y) 

Using Holder inequality, we have 
/ %,...,p}0u(y)c£u 



/ %,...,p}0 u (y)du. 



< 
< 



• tip- 1 



/ |<9{i,.... p }</>u(y)| 7 du 

J[s,t] 



[s,t]| 7 sup |d{i,..., p }0u(y)| 7 

te[-T,T]p 



so that 



\\<h 



17 < 



sup 9{i,..., P }0t(y) 

te[-T,T]p 



n 1* 



In general, if s < t, let I be the set of indices such that s.; < ti. We show similarly 
that 



ll^[.,t]ll^< 



sup dj4n(y) 

te[-T,T]p 



□ 



3.4. Holder- regularity when a — 2. In this section, we prove Proposition 12.31 where 
Holder-regularity is stated for the limit ^(Mt) (resp. J n (A*t)) when G has finite 
variance (resp. finite moments of any order). In particular, we set a = 2 and we 
assume (3 £ (d, 2d) . Holder-regularity is proven using again moment estimates for 
increments. However since we have not found in the literature Holder-regularity 
result relying on generalized increments in dimension p > 1, we use standard incre- 
ments Z2(p.t) — -^(Ms) (resp. J a {pt) — Ja(Ps)) that will be controled thank to the 
following elementary observation: suppose the family (</>t)teRp satisfies condition 
(P 7 ), then for all T > 0, for all [s, t] C [-T, T] p , 

ll^-^tll^Cr^Hs-t]]^. (3.10) 

To see this, rewrite a standard increment as a sum of d block increments along axis. 
More precisely, for s,t £ R p , define u* G W by = tj for j < i and u l j — Sj for 
j > i so that, for each 0<i<p — 1, u l and u l+1 only differ from at most one 
coordinate. For 7 > 1, we have 

1=1 
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so that condition (P 7 ) implies 



Us 



Mi; < Ctv 1 ' 1 J2\u-s z \<c 



TP 



Remark 3.6. The trick consisting in controlling the moment of standard increments 
by generalized increments in (13.101) cannot be used in Section 13.31 Indeed, similar 
bounds as in (|3.8[) but for standard increments M p (fi t ) — M p (ji s ) combined with 
(|B~TU)) would yield 



m i \M p { i xt)-M p { i x s )\^}<CC T p y 



n{p) a 



l/ct 



tiii (1+ ^ /d) . 



instead of 



which is far from a bound in ||s — t\\ s with 5 > p r equire d in 



the criterion for tightness relying on standard increments (see Kunita ( 1990l Th. 
1.4.1)). The trick works when a = 2 because in this case, the limit is Gaussian 
and we can artificially increase the exponent in the bound thank to the relation 
between moments and variance, see (13. lip . 



Proof of Proposition [273l 

First point: We assume that G has a finite variance, h £ i 1 (K d ) n L 2 (M. d ) and 
(P 7 ) holds for 7 = 1,2. Using Lemma [A. l\ we have, for n{dy) = 4>(y)dy, 



Var(Z 2 ( M )) = 

< 



\p[ T x.rh] | 2 , ; o drdx 



r l+0 

M|2 A „2<^||^||2M ^ ,||2^„c ^/ o_-l-^ (^r 



t\\h\\iAr^U\\i\\h\\i)a a C 2 r- 

2—0/d 




mum fwmwH 



/3 -d 



|2-/3/d 



l\\h\\i 

2 



1-/3/. 



(2d-p)(fi-d) 

for some finite constant C. Next, the properties (Pi) and (P 2 ) together with (|3 . 10[) 
entail that for all s, t 6 W: 



Var(Z 2 (^ s -im)) < C\\s-t 



\3-/3/d 



We derive easily the moment of order 2n of the Gaussian random variable Z2 {p s 



(^^_lVar(Z 2 ( Ms -, t) y 



(3.11) 



C- 



(2n-l)! 



|(3-/3/d)n 



(n-l)!2™- in 

Using a standard regularity criterion for random fields (see e.g. Kunita ( 1990l . Th. 
1.4.1)), Z2(nt) has, almost surely, Holder-continuous path for all index striclty less 
than 3d ^ 2< /"^" • Letting n go to +00, we obtain that Z^pt) is, almost surely, 
Holder-continuous for all index striclty less than 3 . 
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Second point: With p = p a — pt, Proposition 13.11 entails 



Ck{M p (fj, s - /i t )) 



k(kd-/3) 



< c- 



n(p) h 



y tl! fc 



till 



(3.12) 



Now observe that in the intermediate regime, \(p)p@ /n(p) k remains bounded since 
it converges to a as p — > 0. Furthermore, using properties (Pi) and (Pk) together 
with (|3~T0l) in ([3T2]l . we deduce 



c k (M p (p s - p t )) <C||s-t|| 



l+k-P/d 



Next, recall that the moments of a random variable X are expressed in terms of its 
cumulants by the so-called complete Bell polynomials, i.e. ELY fc ] = Bk{c\(X), . . . , Ck(X)) 
with 

-B fe (ci, . . . , C fc ) = ^ Kn(h, in)c l i ...C% 

i\ +2%2 — —k 

where i\, . . . , iy. are non-negative integers and K^^ii, . . . , if.) are coefficients whose 
explicit (involved) form is not required in our argument. Since we have 

Cl {M p (p s - /i t )) 11 • . .c fc (Mp( Ms - p t )Y* = 

when i\ ^ 0, we can assume, without loss of generality, that i\ = 0. For all s, t £ M. d 
such that ||s — tHoc < 1, we have 

c 2 (M p (p s -p t )) i °...c k (M p (p 3 -p t )) i >° < C||s-t||F- ( ' +1 -™ 

< C||s - t||^H 1 -' 3 / d )( fc /2) 

where we use ^4=2 ^ = ^ anc ^ 2z=2 *i — ^/2 together with ||s — t||oo < 1. We 
deduce 



E 



(Mp(p s - p t ) 



<C\\s- 



\(3-0/d)k/2 



and using Corollary! 

E[(J a ( Ps -pt)) k ] <C\\s-t\\g-W^. 

From the standard regularity criterion for random fields (see e.g. Kunita I ( 199C)l 
Th. 1.4.1)), J a (pt) has, almost surely, Holder-continuous path for all index striclty 
less than (3 - /3/d)/2 -p/k. □ 

3.5. Convergence in distributions space. The results on convergence in distribu- 
tions space are based on Lemma 13.71 below relating the generalized random field 



(m p (p)) 



and the parametric random field 



(M p (t)) 



defined by 



M p (t) = 



M p {d y )dy 



t € 



= M p ((h), 

where pt(dy) = sign(ti) • • • sign(td)l[o,t] (y)dy is the parametric family of signed 
uniform measures given in Example 12.11 Recall that this family satisfies the prop- 
erty (P 7 ) for all 7 > 1 so that Theorem 12.21 indeed holds. We also define the 
continuous random fields on M. d which are the possible limits in C(M. d ) under, re- 
spectively, the large ball regime and the intermediate ball regime: 

Z a (t) = Z a (p t ), J a (t) = Ja(pt), t 6 R d . 
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Lemma 3.7. For all <p e V(R d ), 

M p (0(y)dy) = (-l) d I M p (y) ^ (y)rfy. 
Proof: This follows from successive integration by parts: 

= f M p (6 y )ct>(y)dy 

= L (f • • • f ^ -.«<>**••■ **) rfy 

= (-W M p (y) g^- (y)dy. 

□ 

We prove now Theorem 12.41 The proof relies on Lemma 13.71 and on Theorem 
12.21 stating that M p converge in C(M. d ) to Z a (resp. J a ) in the large ball (resp. 
intermediate ball) regime. 

Proof of Theorem 12.41 

First point: The random field M p is a bounded linear operator on V{R d ). Indeed, 
for cj) S T>(R d ) whose support is included in [— T, T] d , Lemma \3 . 71 entails : 

Mp(<Xy)dy) < (2T) d ( sup \M P \}( sup — J . 

\[-T,T] d ) \[-T,T] d oy\--- oy d J 

This proves the continuity of the linear application on T)(R d ) and hence that M p 
can be seen as a random distribution. 

Second and third point: Let / : C(R d ) — > V(R d ) be the canonical injection 
given by 

/(/):0h> f /(y)0(y)dy, f eC(R d ), 
and define the differential operator D : V'(R d ) -> 2?'(R d ) by 

D(s):cl ) ^(-l) d s( —^- Y seV'(R d ). 

With these notations, Lemma T3.7I states that M p = (D o I)(M p ). Since, the oper- 
ators D and I are continuous, and M p converges in C(R d ) as p — >• to Z a (resp. 
J a ) in the large ball regime (resp. in the intermediate ball regime), the continuous 
mapping theorem implies that M p weakly converges in T>'{R d ) to {DoI)(Z a ) (resp. 
to (D o I)(J a )). Finally, since weak convergence in V(R d ) implies fdd convergence 

on V{R d ), we have (D o I){Z a ) f = Z a and (D o I)(J a ) f = J a . This shows that Z a 
and J a have modifications that are continuous on T>(R d ). □ 
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Appendix A. Technical results 



Lemma A.l. Let 7 > 1 and n(dy) = (j)(y)dy with (f> G L 1 (R d ) n i 7 (R <i ) and 
ft, G L 1 (R <i ) n L 7 (R d ). Then 



W^h]Vdx < (r d ||</>||7||ft|| 7 ) A (r 7d |H| 7 ||ft|| 7 ). 



(A.l) 



Proof: Since ft, G L 7 (R d ) and /i(dy) = 0(y)dy with G L 1 (R d ), the Holder 
inequality entails: 



y - x 



< ll^lir 1 

= ^ii^iir 1 



<Ky)dy 
y x 



dx 



|0(y)|dydx 



|ft(y)| 7 |<Xn, + :z;)|dydx 



\J\\hU- 



(A.2) 



On the other hand, still using Holder inequality but with <fi G L 7 (R d ) and ft, G 
i" 1 ^), we have: 



|/i[r Xir /i]| 7 dx 



= r 7d 

< r 7d 

< r 7rf|U||7 



7-1 

IMy)My' 



rfx 



|0(ry + x)| 7 |/ l (y)|rfyrf> 



7 

7ll'»lll' 



The bounds ljA.2|) and (|A.3|) together entail (|A.1|) . 



(A.3) 
□ 



The following r esult proves the co ntinuity required to apply Lemmas 2 and 3 instead 
of Lemma 6 in lKaj et al. I ( 2007 ) in the modification in Section I37T1 



Lemma A.2. Suppose that the fading function ft satisfies ([A3]). For /i(dy) G 
L 1 (M d )CiL a (M d ), the application r \-t J Rd ^ G (fi[T x , r h\) dx is continuous on (0, +00). 
The same holds true for J Rd f Q (ti[T x , r h]) dx with ^ a {6) = -cr Q |6»| Q (l+ibe(0) tan(?ra/2)) 
and a, b given in flAil ). 



(A.4) 



Proof: Observe first that, for r n — > r$ > 0, n — > 00, we have: 
lim n[T Xt r n h] — /i[T x .,. Q /i] dx-a.e. 



This is a standard application of Lebesgue's convergence theorem. Indeed, since ft 
is almost-everywhere continous, we have T x , rn ft.(y) — > T x ro /i(y), n — > +00. From 
the definition of ft*, the convergence is bounded for all n > 1: 

Kr n h(y)<f>(y)\ < r x , R h*(y)\cf>(y)\, x,yeR d 

with R = sup{r„; n > 1} and the bound is dx-integrable since Lemma I A . 1 1 applied 
to ft* G n L a (R d ) rewrites 



Tx, R h*(y)\<j>(y)\dy 



dx < 
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A second application of Lebesgue's convergence theorem yields 

lim / ^ G (fi(T x _ rn h))dx = / i& G (ji(T x , rn h))dx. 
n-^+oc j Rd j Rd 

The convergence (|A.4j) together with the continuity of ^>g imply indeed the point- 
wise convergence 

lim ^ G (n(T x .r n h)) = * G (^(r x . ro /i)) dx-a.e. 

n— >+oo 

The convergence is bounded by C\/j,(T Xj Rh*)\ a since \E'g( m ) < C| u | a j which is in- 
tegrate by Lemma ED applied to h* € L 1 (M <i ) n L a {R d ). A similar proof holds 
for 

lim / V a (jj,(T Vn h))dx= / y a (n(T Xtrn h))dx. 

n^oo J Rd J Rd 

□ 

The following result estimates the truncated moments of F. In particular, the 
condition > d ensures that F has a finite moment of order d. 

Lemma A. 3. For 8 > 0, when u — > +oo, we have: 

f Cst if<5</3 
/ r 5 F(dr)~^ /3C/3 In u if 8 = (A.5) 
J ° 1 T^Cpu s -P H8>0 

and for < 5 < 0, when u — > +oo, we have: 

Moreover when S > 0, we have the global bound 

r 5 F(dr) < Cu 4- ". (A.7) 
Proof: Let i? be a random variable with distribution F. We have 

/ r s F(dr) = E [ii 5 l fl < u ] = / / l {R < u} l {t < R s } dtdP 
Jo if!JI+ 

r 6 

"(t 1 ^ <R< u)dt = / P(t 1 /' 5 < i? < u)rft 



(A. 



= 5 P(s<R< w)s* -1 ds 
./o 

= 5 / P(i? > s)sp- x da -8 1 P{R> uJs'-'-ds. 
But using condition QA2D, we have: 



.Cst if<5</3 
/ P(E > s)s s ~ 1 ds ~ ^ C/jlnu if 8 = 
Jo 1 ^u*-" if 8 > 
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and 5 J? P(R > u)s s 1 ds ~ Cpu 5 @ from which (IA.5I) easily derives. Next, 

/ r s F(dr) = E[R S 1 R > U ] = / / l {R > u} l {t < RS} dtdP 
Ju JnJR+ 

F{R > max(t 1/s ,u))dt = 8 1 F(R > max(s, u^s^ds 



= S 



/ P(R > u)s 5 - 1 ds + 5 P{R > s)s 5 - 1 ds 

JO Ju 

which is (|X6l) . Finally, since P(R > s) < 1, (|X8|) entails J^r 5 F(dr) = 0{u s ) so 
that together with (IA.5I) , it is easy to derive (|A.7|) . □ 

Lemma A.4. Let 7 > /3/d and fi(dy) = cj>(y)dy with cp e L 1 (M d ) n L 7 (R d ) and 
/i e L 1 (R d ) n L 7 (R d ). Then, for any p > 0, 

xr+ \l a ~P)\p-a) 
where M is a constant depending only on F. 

Proof: Using Lemma \A.1\ we have: 

M T x,r-fr]P 'dyiF p (dr) 

R+ 

(A d ^ll7l^ll^)A^V d ||0|| 7 ||/,||?) F(dr) 

j-c/p p+OO 

= ll^ll 7 ,^!!^ / r^F{dr) + UWIWHY / r d F(dr) 

JO Jc/p 

with c = ((ll^llill^lUj/Gl^lhll/illi)) '-'i v > i ■ Finally using the bound on the truncated 
moments of F in Lemma [A.3I with S :— -fd > (3 in (|A.5I) and S :— d in (IA.6|) . we 
derive 

/ | M [t x ,^]| W p (dr) < MpP [^^|0|| 7 ||/,|| 7 C 7 ^ + -l-W^JWhW^ 

where M € (0, +00) depends only on _F. The result is obtained, after cancellation, 
by replacing c by its definition. □ 

The next result collects explicit formulas for the characteristic function and the 
cumulants of the rescaled and centered random vari a ble M (p). It is based on 
standard results for Poisson integrals, see iKallenberg ] (l2002h . 

Lemma A. 5. (1) The characteristic function of M p (p) writes: 

Vm (JP) = ex P ( / *G (nipy^piTx.rh}) \{p)d*F p {dr) 

where * G (n) = J R {e lum - 1 - ium)G{dm). 
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(2) Suppose G has a finite moment of order k > 1, fJ-(dy) = 4>(y)dy with 
(f> G L 1 (R rf ) n L fe (R d ) and /i G L 1 (E d ) n L k (R d ). Then M p (^t) has a finite 
moment of order k and its k first cumulants are given by: 

ci(M p Qj)) = 0, 

q(M») = ^4 ( I ml G ( dm )) ( I W*,rh]) l dx F p (dr)) ,2<l<k. 
n(p) \Jm J VjR d xR+ / 

Note that the finiteness of J R d xR + iMlTx.r'i]! dx.F p (dr) comes from Lemma [A. II 
when n 6 L 1 ^) H L k (R d ) and ft G L x (K d ) n L fe (R d ). 

References 

P. J. Bickel, M. J. Wichura. Convergence criteria for multiparameter stochastic 
processes ans some applications. Ann. Math. Stat, vol. 45, no. 5, pp. 1656-1670, 
1971. 

H. Bierme, A. Estrade, I. Kaj. Self-similar random fields and reseated random balls 
models. To appear in J. Theoret. Probab., 2010. 

J.-C. Breton, C. Dombry. Rescaled weighted random balls models and stable self- 
similar random fields. Stoch. Proc. Appl., vol. 119, pp. 3633-3652, 2009. 

J. Deshayes, D. Picard. Principe d'invariance sur le processus de vraisemblance. 
Ann. I.H.P., Sect. B, vol. 20, no. 1, pp. 1-20, 1984 (in french). 

W. Feller. An Introduction to Probability Theory and its Applications. Vol. 2. Wiley. 
1966. 

I. Kaj, L. Leskela, I. Norros, V. Schmidt. Scaling limits for random fields with 
long-range dependence. Ann. Probab., vol. 35, no. 2, pp. 528-550, 2007. 

I. Kaj, M. Taqqu. Convergence to fractional Brownian motion and to the telecom 
process: the integral representation approch. In Brazilian Probability School, 10th 
anniversary volume, Eds. V. Sidoravicius, M.E. Vares, Birkhauser, 2007. 

0. Kallenberg. Foundations of Modern Probability. 2nd Ed. Springer New- York, 
2002. 

1. Kaj. Aspects of wireless network modeling based on Poisson point processes. 
Fields Institute Workshop on Applied Probability, Carleton University, available 



at http : //www. math. uu. se/~ikaj /preprint s/wirelessnotes . pdf , 2006. 



H. Kunita. Stochastic flows and stochastic differential equations. Cambridge Stud- 
ies in Advanced Mathematics, vol. 24. Cambridge University Press, Cambridge, 
1990. 

T. Mikosch, S. Resnick, H. Rootzen, A. Stegeman. Is network traffic approximated 

by stable Levy motion of fractional Brownian motion. Ann. App. Probab., vol. 

12, no. 1, pp. 23-68, 2002. 
B. von Bahr, C. G. Esseen. Inequalities for the r-th absolute moment of a sum of 

random variables, 1 < r < 2. Ann. Math. Statist, vol. 36, pp. 299-303, 1965. 
X. Yang, A. P. Petropulu. Co-Channel interference modeling in a Poisson field of 

interferers in wireless communications. IEEE transactions on signal processing, 

vol. 51, no. 1, pp. 64-76, 2003. 



